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Abstract 

Integrate & fire (IF) neurons have found widespread applications in computational neuroscience. 
Particularly important are stochastic versions of these models where the driving consists of a synaptic 
input modeled as white Gaussian noise with mean ^ and noise intensity D. Different IF models have 
been proposed, the firing statistics of which depends nontrivially on the input parameters ^ and D. In 
order to compare these models among each other, one must first specify the correspondence between 
their parameters. This can be done by determining which set of parameters [ji, D) of each model is 
associated to a given set of basic firing statistics as, for instance, the firing rate and the coefficient of 
variation (CV) of the interspike interval (ISI). However, it is not clear a priori whether for a given 
firing rate and CV there is only one unique choice of input parameters for each model. Here we 
review the dependence of rate and CV on input parameters for the perfect, leaky, and quadratic IF 
neuron models and show analytically that indeed in these three models the firing rate and the CV 
uniquely determine the input parameters. 



1 Introduction 

Stochastic integrate & fire (IF) neurons constitute an important tool in theoretical neuroscience, hav- 
ing been used to address a number of relevant biological problems. For instance, different variants of 
these models have been empl oyed in the debate on the high variability of the inters pike interval (ISI) 
observed for cortical neurons (j'Softky and Kochl . 11993 IGutkin and Ermentroutl [T998f ). O ther problems 
in which stochastic IF models have been applied include the response to fast signals (IBrunel et all . 
[20011 ILindner and Schimanskv-Geieii . [200l IFourcaud-Trocme et al.l . [2001 INaundorf et al.U2005l ). asyn- 
chronous spiking in recurrent net works (IBrunel , 2000 1) , and oscillations of f iring ac tivity in systems with 



spatially correlated noisy driving (jPoiron et al.l . [200 3. 20041: [Lmdner et "all . 12005 b!). 

IF neurons can be classified according to the nonlinearities that govern their subthreshold dynamics. 
Three simple and important variants are the perfect (PIF), leaky (LIF), and quadratic (QIF) models, in 
which the subthreshold dynamics is described by a constant, a linear, and a quadratic voltage dependence, 
respectively. Noisy inputs with different degrees of biological realism have been considered for these 
models. One simple choice is a white Gaussi an inpu t curr e nt, corresponding to the so-calle d diffusion 
approximation of synaptic spike train input ([Holdenl . I1976t IRicciardil . I1977t iTuckwelll . 119891 ). The PIF 
wit h white noise drive (also referred to as the random walk model of neural firing) was first considered 
by IGers tein and Mand elbrotI (| 19641 ). The LIF with a white Gaussia n input curren t has been studie d 
by (IJoh anncsma. 1968) and a fterwards by many other authors (see IHolden[ (|1976D : IRicciardi[ (|1977f ): 



ITuckwelll (il989i):iBurkitti (120061) and references therein); it is also referred to as Or nstein-Uhlenbeck neuron 
( Lanskv and Rosparsl. I1995l)~^ W hite-noise driving in the QIF was studied by IGutkin and Ermentroutl 



19981) and lLindner et al.ll[2003f) . 



The firing statistics of the various IF models may depend sensitively on the specific nonlinearity in the 
respective model , and so it is not clear at the first glance which model is capable of reproducing which 
features of the firing statistics of real neurons. For instance, it has been shown that LIF and Q IF display 
rather different phase shifts if driven by a periodic stimulus (jFourcaud-Trocme et al.l 120031 ) . Further, 
the LIF with perio dic stimulation can transmit signals o f very high frequencies if they are encoded i n 
the noise intensity ([Lindner and Schimanskv-Geieii 120011 ). while the QIF cannot ([Naundorf et al.l . 120051 ). 
Last b ut not least, when the s t rength of the input n oise is varied, the LIF can show coherence resonance 
(CR) ([Pakdaman et all . [200l ILindner et all . [2002 ) whereas the PIF and the QIF do not display CR 
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(|Lindner et al.l . I2003D (see also below) . A thorough comparison of different stochastic IF neurons and the 
roles of their respective nonlinearities constitutes therefore an interesting and largely open problem. 

If one wants to compare two specific IF models one should tune their input parameters (mean and 
intensity of the input fluctuations) such that the basic firing statistics is the same. A simple choice for 
the basic firing statistics is the firing rate, quantifying the intensity of the spike train, and the interspike 
interval's coefficient of variation (CV), quantifying the irregularity of the spike train. Setting both models 
in this way in the same firing regime (e.g. low spike rate and high ISI variability), one can then ask for 
higher statistics (e.g. the power spectrum of the spike train) or for their response characteristics (e.g. to 
weak periodic stimulation or step currents). The idea for such a tuning tacitly assumes that there is for 
each of the models at most one input parameter set yielding a desired firing statistics, for instance, rate 
and CV. At a closer look, however, this is not evident at all: are the input parameters for a white-noise 
driven IF model indeed uniquely determined if we prescribe certain values of the rate and the CV? This 
is the question that we address in this paper and we will answer it for the three IF models mentioned 
above, namely, PIF, LIF, and QIF neurons. 

The question discussed here is also related to the problem of par ameter estimation for a n IF rn odel from 
experimental data, which has been subject of several studies (see ILanskv and DitlevsenI (j2008D and ref- 
erences therein). In one approach, t he model pa rameters are inferred from subthreshold membrane mea- 
surements (for a recent reference see lBadel et al. (2008) ). In another approach, model parameters a.re es- 
timated using solely the ISI statistics (ITuckw'eiTand Richte?, TO TSlMInoue et all [19951: IRauch et a l.'. '20031; 
Camera et all . 120041: IShinomoto et all . 119991 : IDitlevs en and Lan skvl . 120051 l2007l:lMullownev anT lvcngar. 



2008( ). The latter approach is of importance, since often subthreshold data are not available. Most of 



these studies consider the leaky IF model. To th e best of our knowle dge, the problem of the uniqueness 
of input parameters has been addressed only bv IKostal et al.l (|2007l '). where it is mentioned without a 
proof that uniqueness of the input parameters given fixed rate and CV holds for the LIF. 

In this paper, we show analytically that rate and CV uniquely determine the input parameters for 
the three models. We first introduce the models and the statistics studied here (sec. [2]). In the main 
part of the paper (sees. [3][5]), we briefly review the rate and CV as functions of the input parameters and 
then prove the uniqueness of the relation between these parameters and the former statistics. We give 
an outlook to generalizations of the considered problem in sec. |6l 



2 Integrate & fire neuron models 

2.1 Definition of the models and relation to the first passage time problem 

IF models consist of two ingredients: (i) a one dimensional stochastic ordinary differential equation 
describing the subthreshold behavior of the membrane potential 1^ as a function of time t and (ii) a 
fire-and- r eset r ule. The equation for the membrane potential has the form of a current-balance equation 
(|BurkittJ . [20?lh : 

= L,.UV) + lUt) + lUt), (1) 

where Cm is the membrane capacitance, /ay„(i) and I^„t{t) denote the synaptic and injected current, 
respectively. Here we will consider the case where I(t) = Isy„it) + /oxt(i) is a white Gaussian noise current 
with a constant mean value (/) and a correlation function {{I{t) — {I)){I{t') — (/))) = 2DjS{t — t'). The 
function ImodeiiV) stands for a model-specific current — for instance, a passive leak of the membrane 
would be described by setting Imodei = — Vg] / Rm , where Vq and Rm denote the leak reversal potential 
and the passive membrane resistance constant, respectively. 
The fire-and-reset rule can be expressed as 

V{t) = Vth spike at time t and V ^Vr, (2) 

i.e., whenever the membrane potential reaches a threshold value Vth the neuron fires a spike and there is 
a reset of its membrane potential to a value Vr- 

It is convenient to make the following changes of variables: 

i^^, (3) 



Vth - Vr T„ 

where t,„ = CmRm is the membrane time constant and the new variables, v and t, are dimensionless. 
This procedure corresponds to measuring the voltage in units of the difference between threshold and 
reset (with Vr as the reference voltage) and time in units of the membrane time constant. 
Defining 

/„o.. = 777^^^ [L..UiVth - Vr)v + Vr) ~ /.„„..,(K)] , (4) 
[Vth — Vr) 
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and 



we can recast Eq. ([T]) into the form: 



JTT^^J{I)+I....AVr)], (5) 
[Vth — Vr) 

D R'^ 

D= J (6) 



v = U.M+^^ + V2Dat), (7) 



which is the equation that we will work with in the remainder of this worlQ- The parameters ^ and D 
are called input parameters, i.e. they represent the mean and the intensity of the fluctuating input in 
our nondimensional model. Note that they are linearly related to the physiological input parameters (/) 
and Dj via Eq. ([5]) and Eq. ([6]). The function ^{t) is a zero-mean white Gaussian noise with {£,{t)(,{t')) = 
6(t — t'). In the nondimensional formulation Eq. ([7]), the reset and threshold values. Vj. and Wj^, are 
given by zero and one, respectively. For better applicability of our results, we will keep Vr and vth in all 
resulting formulas and equations. 

If one neglects any voltage dependence of the right-hand side in Eq. ([7]), one deals with a perfect 
integrator and the respective IF model is called a perfect integrate & fire (PIF) neuron: 

/pip = 0, /i > 0, Vth = 1, Vr = 0. (8) 

If a leak current is taken into account, f{v) becomes a linear function (by construction, any additive 
constant is lumped into the mean input /i). In this case, we deal with a leaky integrate & fire (LIF) 
neuron which is the most often studied and used IF model to date: 

/lif = -W, /iG (-00,00), Vth^l, Vr^Q. (9) 

Another IF model, the quadratic integrate & fire (QIF) model also fits into the framework of Eq. ([7]) 
although in its derivation v is not a voltage but a more abstract variable. Higher nonlinearities in the 
current-balance Eq. ([T]) come along by voltage-dependent conductances and additional variables like gat- 
ing variables (e.g. in the classical Hodgkin-Huxley model) or effective recovery variables (e.g. in the 
Morris-Lecar model). A wide class of neurons (type I neurons) described by su ch multidimensional dy- 
nami cal systems shows a behavior associated with a saddle-node bifurcation (jRinzel and Ermentroutl 
119891) in these multidimensional models. For these neurons, the dynamics close to the bifurcation is gov- 
erned by one slow variable only which is the variable v in our general IF model Eq. ([7]). The nonlinearity 
corresponding to the normal form of a saddle-node bifurcation is a quadratic function. Furthermore, for 
this specific dynamics the behavior is governed by the slow motion close to w = and threshold as well 
as reset values do not matter if taken far away from the origin. For simplicity, they are commonly taken 
at infinity. Thus the dynamics of a QIF model is described by the following setup: 

/qif = W^, /^e (-00,00). Vth^QO, Vr = -00. (10) 

It is useful to interpret Eq. ([7]) as describing a Brownian particle of position v undergoing overdamped 
motion in a potential [/modoi such that 

atVmodcl 

-—; =/„„de,+M. (11) 

av 



In this analogy, the ISI of the respective neuron model turns into the first-passage time of the Brownian 
particle starting at the reset point towards the threshold point. Depending on the model and, in particular, 
on the value of /i, the passage can occur already without noise {tonic firing regime) or must be assisted by 
fluctuations (noise-induced firing regime). The tonic regime can be most easily illustrated in case of the 
PIF where the particle just slides down an inclined plane from reset to threshold (see Fig. [T^) , whereas 
noise is needed to reach the threshold whenever there is a barrier present between reset and threshold 
(QIF for n <0, see Fig. (T};) or right at the threshold (LIF, n < 1, see Figs.[He) and (f)). Note that the 
parameter fi has different meaning in the three models. In the PIF it attains only positive values and sets 
merely the time scale of the system. In the QIF it is a bifurcation parameter: at negative /i the potential 
attains one minimum whereas for positive fi the potential is a nonlinear but monotonic function. In the 
LIF, the bifurcation from tonic to noise- induced firing takes place at /i = 1. As we will see, for the firing 
statistics of the LIF it is furthermore useful to distinguish the case where /i < 0: here large values of the 
CV can be easily interpreted in light of the specific properties of the potential (see below). 

^Note that, by construction, /modoi(ffl) = 0. Another choice for v is possible that also leads to Eq. (0 in which v 
denotes the deviation of the membrane potential from the leak reversal potential Vq instead of the scaled deviation of the 
membrane potential from the reset potential Vr as in our case. 
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Figure 1: Potentials for the different models (of. Eq. (fTT|) ). The PIF (a) is only defined in the tonic firing 
regime (/i > 0) while the QIF ((b) and (c)) can be in the tonic (QIF I in (b)) or in the noise-induced 
(QIF II in (c)) firing regimes. The LIF displays tonic firing (LIF I in (d)) for ^ > 1 and two different 
noise-induced firing regimes ((e) and (f)) for which we distinguish the cases where the potential minimum 
lays between and Vth (e) and where the minimum is to the left of the reset point (f). In the latter 
case, the potential picture allows for a simple interpretation of the high variability {CV > 1) under 
moderate noise intensity: short ISIs occur when the particle moves directly from v^. to Vth, while long 
ISIs correspond to realizations in which the particle first performs an excursion to the potential minimum 
and then hits threshold. 



2.2 Measures 

The spike train is defined as a sum of delta functions at the spiking times (jGerstner and Kistlerl . [200^ . 
i.e., the time instants when the voltage reaches the threshold and the fire-and-reset rule is applied (cf. 
Fig. ID): 

j/(t)=^<5(t-t,). (12) 

j 

In Eq. (fT2|) , tj stands for the instant when the j-th spike is triggered. Fig. [2] depicts the time evolution 
of the subthreshold voltage as described by one of the models we address and the corresponding spike 
train. The time intervals Tj = tj — tj-i between two immediately subsequent spikes are the ISIs. 

spike train 



'O <2 'S <4 

time 

Figure 2: Subthreshold voltage dynamics and corresponding spike train as from a simulation of Eq. ([7]) 
and Eq. ([8]) (PIE) with parameters /i = 0.9 and D ~ 0.006. Note the high degree of regularity of the 
spike train, characteristic of the PIF model under strong mean input with low noise intensity. 



The spike trains considered here are stationary stochastic point processes. The firing rate r of such a 
process can be defined as the inverse mean ISI: 

The coefficient of variation {CV) of the ISI is defined as: 



cv^^. (») 



4 



where (AT^) = (T^) — (T)^ is the variance of the ISI distribution. The CV can be regarded as the relative 
standard deviation of the ISI. For comparison, a perfectly periodic spike train would have zero CV while 
a Poissonian spike train possesses a CV of one. 

As regards experimentally measured values of the rate and the CV, wc note that all nondimensional 
values of rates discussed in the following translate into the former by r^ea; = r/r. For example, for a 
membrane time constant r = RmCm — 10ms nondimensional IF rates of r = 0.01,0.1, and 1 correspond 
to firing rates of rreai = 1, 10, and 100 Hz. Note that the CV (being a relative standard deviation) does 
not change upon a time transformation, i.e. CVreai = CV. 



2.3 General form of the differential equations governing the contour lines 

Analytical formulas for the mo ments (T") of the first passage time from x- to x+ in an arbitrary 
potential U{x) were derived bv IPontrvagin et al.l (jl933l ). Simplifications of these quadrature formulas 
as well as sur n formu l as for specific cases ha,ve been put forward by many authors (for a selection, see , 
for instan ce. (IHoldenl. \WM IRicciard i [19771 IRi cciardi and SaceTdotl [19791 ITuckwelll . [l98l IColet et al.l . 
[1989: Bul sara et all . 119941 . 11996: Lin dner et al.l . 2002, 2003)). The first two moments determine the rate 
and CV, according to Eqs. ([T3|) and ([T4|) . For completeness, we write the expressions for these moments 



here: 



and 



(T) = 1 r da;e^(-)/^ P dye-^^y^^ (15) 

^ J X- J — OO 

^ Jx— J — oo Jy J—oo 

In this paper we will study the rate and CV of the three models as functions of the input parameters, 
fi and D. In particular, we are interested in the curves for which F{D, ^) — const (F denotes either r or 
CV), i.e., the contour lines of the surfaces F{D,fj,) over the {D,^) parameter plane. The contour lines 
can be obtained if one observes that dF = along them. Since dF = {dF/dD)dD + {dF/dn)dfj,, one 
concludes that the following differential equations for functions firiD) or -Df(m) parametrize the contour 
lines: 

di,F _ dF/dD ^^^^ 



dD dF/dfi ' 

dPp _ dF/dn 
~~dir ^ ^ dF/dD' 



(18) 



where F e {r,CV}, provided that dF/dfi ^ and dF/dD ^ 0, respectively. We note that these 
conditions are not necessarily satisfied in the whole (D,^) parameter space of the models we address. 
For instance, for the PIF we have in fact dr/dD = for all valid pair {D,fi). However, for the three 
models studied here, at any point {D, ji) of parameter space at least one of these conditions is satisfied. 

If for any pair {r,CV) the respective contour lines ^r{D) and ^cv{D) intersect at most once, then 
rate and CV determine uniquely the parameters of the respective IF model. In the following sections we 
will show that this is indeed the case for the PIF, LIF, and QIF. 



3 Perfect integrate &; fire neuron 

The mean and variance of the ISI are given by (jHolded [19761 ITuckwell [19881: [BuTsara et al.l . [1991 : IBurkitE 
[2fM) : 

(T) = (Ar^) = '''("V"'-). (19) 

We stress that ^ > for the PIF; otherwise all moments of the ISI diverge. For this model the expressions 
for rate and CV are quite simple: 

r=—^, CV^= , (20) 

Vth-Vr fl{Vth-Vrj 

Moreover, the contour lines for the rate and the CV can be explicitely calculated (without resorting to 
the differential equations Eq. ([TTjl and Eq. ([T8|l ): 



2D 

firoiD) = ro{vth - Vr), fJ.cVo{D) 7 T77f72- (21) 
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We briefly review the behavior of rate and CV as functions of /i and D and then show that rate and CV 
uniquely fix the system's parameters. 

3.1 Rate and CV and their contour lines in the {D,fi) plane for the PIF 

The rate and CV are shown in Figs. [Sja) and[3l[b) as functions of the parameters fi and D. The rate is a 
hnear function of /i and, remarkably, does not display any dependence on D. This is a unique property of 
the PIF model. The CV depends linearly on y^D/ fi, and can therefore attain values in the whole range 
0<CV <oo. 




Figure 3: Rate (a) and CV (b) as a function of the parameters fi and D for the PIF model. The rate 
is a simple increasing function of the mean input and does not depend on the noise intensity. The CV 
is proportional to y^D/ fi (cf. Eq. (|20p). In (c) we show contour lines corresponding to three different 
values of the rate (dashed lines) and to four different values of the CV (solid lines with symbols). If for 
a given pair of rate and CV there is only one pair of input parameters (/i, D), there should be only one 
intersection point for two specific (dashed and solid) contour lines. The uniqueness of this correspondence 
can be explicitely shown for the PIF (see Eq. ([22]) and Eq. (p3|) ). 

Fig. [3]Jc) shows contour lines for different rates and CVs, which are for both measures just straight 
lines. Generally, the variability of the PIF's spike train increases by decreasing the mean input and 
increasing the noise intensity, which is quite intuitive. 
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3.2 Uniqueness of the model parameters for a given rate and CV for the PIF 

In the fairly simple case of the PIF, Eq. ([20|) can be readily inverted to yield /i and 13 as a function of 
rate and CV: 

= r{vth - Vr), (22) 



D = 

Eqs. dill) and ^ define a mapping (r, CV) 
one and only one pair (D,fi). 



(23) 



{D,^), implying that for any pair (r,CV) there exists 



4 Leaky integrate &; fire neuron 

For the LIF, the mean and variance of the ISIs are (jRicciardil . [19771 : ITuckweI1l . ll988t [Lindner et al.l . l2002f) : 



(T) = / dye^" erfc(2/), 

J a 

(AT^) = 27r / dze^' / d?/e^'erfc^(2/), 



where 



a = — Vth)/V2D and b = {fi — Vr)/V2D. 



(24) 

(25) 
(26) 



From these expressions and the general relations Eqs. pT]) and Eq. (US]), one can derive the differential 
equations that govern the contour lines as follows: 



dD 
F(a,6) 



5- 



' be'^ erfc(fe) — ae°- crfc(a) 
vth — Vr \ e''^ erfc(&) — e"^ erfc(a) 

b — a 



Vth - Vr , 



a{l-¥{a,b)y' + b\ 1 



F(a,5) 



dD 



CV 



dxe^ evic{x)e^ Jl^ dye^ erfc (y) — 2 dze^ dye^ erfc {y)e'' erfc(&) 
JI^ dxe^^CTk{x)e°-^ dyey^evi(?{y) — 2 dze^^ d?/e^^erfc^(?/)e"^erfc(a) 
ducv 



d^jL 



dD 



(27) 
(28) 
(29) 
(30) 



We will fir st recall some properties of rate and CV, mos t but not all of which have been already discussed 
elsewhere (jPakdaman et al.l . 120011 : ILindner et"aI1 . 120021 ) . 

4.1 Rate and CV and their contour lines in the (-D,/u) plane for the LIF 

Rate and CV as functions of /i and D are shown in Fig.Hl As seen in Fig.HJa), the rate is an increasing 
function of /i for fixed D and an increasing function of D for fixed /x. In the zero-noise limit, the rate is zero 
for ^ < Vth and increases monotonically for ^ > vth according to the well-known rate of a deterministic 



LIF model r ^ = ln[(/i — Vr)/{^ 



Vth, 



I IGerstner and Kistleil ([2003 ). 



The behavior of the CV is much richer (see Fig.[4fb)) . In p articular, the LIF model displays coherence 
resonance (CR) (jPakdaman et all . [200lt ILindner et all . [2002: ): for fixed fi < Vth, the CV exhibits a 
minimum at a finite value of D. Coherence resonance thus corresponds to the phenomenon by which 
noise has the counter-intuitive effect of increasing the regularity of the spike train. For the LIF, a 
pronounced CR is observed for a mean input /i close to but smaller than the threshold {fi < Vth)- 

The CV for the LIF can exceed unity (this feature is also displayed by the PIF but not by the QIF). 
Loosely speaking, such a regime corresponds to a firing activity more irregular than in the Poissonian 
regime {CV = 1). This high variability is associated to short ISIs occurring relatively frequently, but 
long ISIs being also likely. When < 0, a simple interpretation can be made in terms of the Brownian 
particle in a parabolic potential. As shown in Fig.[TI^f), in this case, both v^ and Vth arc larger than the 
value of V at which the potential attains its minimum. The short ISIs then correspond to the cases when 
the particle heads directly from Vr to vth, while the long ones correspond to the particle first going to 
the minimum of the potential and then performing its excursion to Vth- 
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Independently of D, if /i ^ — (X3 the firing becomes Poissonian {CV =1). In the opposite Hmit of 
^ oo, the firing is perfectly regular {CV = 0). At least for large noise intensity, the CV exceeds unity, 
as discussed above. Therefore, for fixed and sufficiently large value of the noise intensity we observe a 
maximurr{3 of the CV with respect to /i. This is an interesting feature of the LIF model - not shared by 
PIF or QIF - which to our knowledge has not been described so far. It implies that, given a fixed level of 
the input fluctuations, the degree of irregularity of the spike train is maximal (maximized incoherence) for 
a finite value of the mean input. We note that another type o f maximized incohere nce has been described 
for the LIF with an additional absolute refractory period bv lLindner et al.l (j2002f ). In that case, the CV 
achieves a maximal value for a finite noise intensity when the mean input is fixed. 

As also shown in Fig. [^b), the contour lines for the CV display nonmonotonicities with respect to 
both parameters. The contour lines at which the CV is smaller than 1 display nonmonotonic behavior 
with respect to D, whereas the ones corresponding to CV larger than 1 arc nonmonotonic functions of fi 
(sec, for instance, the contour line CV = 1.1). Fig.HJc) shows contour lines for rate and CV in a range 
of physiological interest. 

We also observe in Fig. lUJd) that the contour lines of rate and CV are very close to each other in 
the region of small D and fi < Vth (see especially contour lines r = 0.01 and CV = 0.95). The reason 
is that in this region we approach the Poissonian firing limit, in which case the dynamics of the model 
depends on the threshold va lue — hence on a — but not on the reset value — hence not on b (see e.g. 
IDitlevsen and Lanskvl (j2005l ) and references therein). For a small but finite value of D and /i < Vth, the 
dependence on the parameter b that distinguishes the contour lines of rate and CV is very weak which is 
the reason why these contour lines will be very close to each other. Assuming an exclusive dependence on 
a, we obtain a universal contour line for all ISI statistics (this holds strictly true only in the Poissonian 
limit) by setting a =const . One such line (with a = —2.1 extracted from the rate curve for r = 0.01 at 
D = 0.001) is shown in Fig. Iljd): indeed, the corresponding contour lines of rate and CV are very close 
to this line. 

Thus, as the firing regime approaches the Poissonian limit, the actual determination of the intersec- 
tions of the contour lines becomes a practically more difficult task. Also it becomes less clear whether 
there is only intersection point or not. In view of this particular (numerical) uncertainty, but also in view 
of the nonmonotonic behavior of the contour lines ficv a-s functions of D (for /i < Vth) and /i (for strong 
noise), it is desirable to gain certainty about whether rate and CV uniquely determine D and /i in the 
LIF model. 

4.2 Uniqueness of the model parameters for a given rate and CV for the LIF 

Our strategy to show that the model parameters are uniquely determined for a given rate and CV 
comprises two steps. First, we demonstrate that each contour line for the rate is unique. Second, we 
prove that the CV is a monotonic function along any rate contour line. The second step can be simplified 
by noting that the CV is the ratio between the square root of the variance = (AT^) and the mean 
(T). Since the mean is invariant in any contour line for the rate, it suffices to show that the cr^ is a 
monotonic function along any such contour line. In other words, it suffices to show that the directional 
derivative of along the tangent of the contour line for the rate is strictly positive 0. 

4.2.1 Uniqueness of contour lines for the rate 

Let us prove that the contour line for a specific value of the rate r{D,fj,) = Tq is one single connected 
curve. The proof comprises three steps: First, for any point (D, ^) G x IR we can locally construct a 
contour line ^ra{D) parametrized by D such that r(D, ^ro{D)) = tq. This is possible locally by virtue 
of the implicit function theorem since, as shown in Appendix iBl dr/d^ > for all {D,fi) G R+ x M. 
Of course, the specific value of the rate, ro, will depend on the point {D,^). Second, we can extend 
{D) to the whole domain D S R+ by connecting neighborhoods of this domain. This could be made 
impossible if firo{D) diverges at finite D. We rule out this possibility by noting that it is not consistent 
with the limit values of the rate at /x = ±oo. Indeed, as we show in Appendix [Bl the limit of the rate is 
zero (for ^ — > — oo) or cx3 (for fi — > oo), values which are not attained by the rate for any positive D and 
finite /i. Hence no contour line starting within the domain can approach the boundary ±oo at finite D 
and thus there cannot be a divergence of the contour line at finite D. Hence firo{D) will describe a single 
connected line for the whole domain D E M+. Third, we show that the graph of firo{D) contains all 
points [D^fi) with r{D,fj,) = ro. In fact, if a point {D* , ^*) not belonging to the graph of jiroiD) exists 
such that r{D*,^*) = ro, then the condition dr/dfi > must necessarily be violated along the vertical 
line D = D* . This completes the proof that the contour lines for rate are single (connected) curves. 

^The existence of this maximum was pointed out to the authors by Tile Schwalger, Max- Planck- Institut fiir Physik 
Komplexer Systeme (Dresden). 

•^The basic idea for this step in the proof is due to Dr. Jochcn Brocker, Max- Planck- Institut fiir Physik Komplexer 
Systeme (Dresden). 
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Figure 4: Rate (a) and CV (b) as a function of the parameters (x and D for the LIF model. In (c) we 
show contour Unes for three values of the rate (dashed lines) and four values of the CV (solid lines with 
symbols). Note the nonmonotonic character of the CV's contour lines fj,cvo{D). The rectangular region 
marked in (c) is magnified in (d) . Also plotted in (d) is one curve (brown thick line) of the type a = const 
[cf. Eq. ([26]) ]. corresponding to a contour line for the ISI statistics in the Poissonian firing limit. Contour 
lines in this limit for the rate (r = 0.01, black dashed line) and CV {CV = 0.95, magenta triangles) are 
very close to this line and, consequently, to each other. 

4.2.2 Proof that cr'^ is a monotonic function along the rate contour Unas 

Along a contour line of the rate, the mean ISI is fixed by definition and thus the CV can only vary due 
to changes in the variance a^. Thus if we show that the variance increases monotonically as we move 
along the contour line in the direction of increasing D, wc will have also shown that the CV increases 
monotonically if wc move along the contour line in this direction. 
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The monotonicity of cr^ along the rate contour Hnes is expressed by 

Vcr^ • Vt > 0, (31) 

where Vt is a vector which is tangent to these contour hnes and V denotes the gradient in {D, fi) space. In 
order to show that Eq. ([3T|) holds, let us first determine vt. Along the rate contour lines, the differential 
equation Eq. (|17p with F = r holds true; its rig ht-h and side is needed for an expression of the tangent 
vector of the rate contour lines appearing in Eq. pTjl : 

'du' 

where eu and are the respective unit vectors. The relation to be shown, Eq. ([3T|) . thus corresponds to 



= + -^T^e^ (32) 



It is much simpler to express the derivatives on the left hand side of Eq. ([33)1 in terms of coordinates 
a and b rather than D and /.i. Using Eqs. ((25|) and ([26|) . we obtain 

dB {Vth - VrJ' 

and 



/OO /"OO \ 

d2/e^'erfc^(?/)+ae"'y dye^' eiic' (y) j (34) 



da^ = Mb- a) ^yg,= g^f^2(y) _ r . (35) 

Vth ~Vr \ Jt J a J 

Inserting Eqs. ((27l [34ll35|) into Eq. (l33ll . and performing straightforward algebra, we write the latter as: 

2^6 - a)^e°^+^%rfc(a)erfc(&) / . eric'jy) ^ . erfc^\ 

(t;t^.-^;.)2(e-^erfc(a)-e^^erfc(6)) U ^ crfc(a) J, ""^^ erfc(6) J ^"""^ 

This inequality holds true for all a < 6 since the two functions erfc(a) and dx erfc^(a:;)/crfc(a) 
(differences of which appear in Eq. (|36p) a re monotonically decreasing functions of a as we prove in the 
Appendix |X1 We have thus proven Eq. (|33|) and hence we have shown that the variance of the ISI always 
increases as we go along the constant rate contour line in the direction of increasing noise intensity. This 
completes our proof of the uniqueness of parameters determined by prescribed values of the rate and the 
CV. 



5 Quadratic integrate & fire neuron 

For the QIF, one has ([Lindner et al.l . [2003D : 



(T) = 



D 



1/3 



/(a), /(a) 



dy e 



ay+y 



(37) 



{AT' 



D 



2/3 



dxe 



a 



3 

D' 



dye-'^y-y 
1/3 

pi. 



dze 



az+z" 



(38) 



(39) 



For this model, the following scaling relations (jLindner et al.l . 120031 ) facilitate the determination of the 
contour lines in parameter space for rate and CV: 



m 

IMI 



(40) 
(41) 
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The scaling relation Eq. (|¥T|) . together with the monotonicity of the CV for /i = ±1 fsee ILindner et ahl 
((20031), implies that (for /i = ±1) a certain value of the CV (say, CVq) determines uniquely the noise 
intensity D which we call D. With this observation, the contour lines for the CV for arbitrary fi can be 
explicitely written: 



D 



2/3 



"■•(^'^RUJ . (42) 

For the rate, D can be regarded as a parameter of the curve ^r{D): from the above definition of D 
and from the scaling relation for the rate Eq. ([40]) we can infer that 

describe all the points on the curve iJ,r{D) which we get by varying D. 

We will now recall some properties of rate and CV which have been already discussed bv lLindner et al.l 
(EOQl. 

5.1 Rate and CV and their contour lines in the (D,fi) plane for the QIF 

Rate and CV as a function of yu and D are shown in Fig. [S] The behavior of the rate is similar to the 
case of the LIE. Here again it is a monotonically increasing function of /i for fixed D and monotonically 
increasing function of D for fixed fi. A noticeable difference arises in the zero-noisc limit and close to 
the bifurcation at /i = 0: the rate for the QIE is also strictly zero if /i < but, differently from what is 
observed for the LIE, it increases proportionally to the square root of /i for small positive fi. 

In clear contrast to the LIE and PIE, the CV for the QIE is bounded in the interval < CV < 1. 
Moreover, here the CV is a strictly monotonic function of both D and /i. Eor fixed D, it decreases with 
increasing /x and, for fixed positive (negative) /i, it increases (decreases) with increasing D. Therefore, 
coherence resonance does not occur for this model. 

Fig-Ht shows some contour lines for rate and CV in the physiologically relevant region of the parameter 
space. We emphasize that for the QIE the contour lines of the CV are monotonic functions of both D 
and fj,. 

As shown in Eig. [5ji, the contour lines of rate and CV are in close vicinity in the region of small D 
and /i < 0, analogously to what is observed for the LIE. In this Poissonian firing li mit, the dyn amics of 
the QIE depends most strongly on the ratio of potential barrier and noise intensity ([Lindner et al.. . ,2003i ) 
which can be interpreted as a rescaled barrier 

Au = (44) 

Eor instance, the firing rate in the Poissonian limit (corresponding to the Kramers rate out of the cubic 
potential we ll) depends on this b arrier exponentially and has only a mild additional dependence on /i via 
a prefactor (| Lindner et al.l . 120031 ) 

r = VH/Tre"'^", Au > 1. (45) 

The line Au =const is indeed close to both the contour lines of rate and of CV in the Poissonian limit (cf. 
brown line in Eig. [SJi). Note that Eq. (j45|) offers another, more accurate approximation, to the contour 
line of the rate in the Poissonian firing regime. This line is, however, very close to the line obtained from 
a constant barrier Au =const. 



5.2 Uniqueness of the model parameters for a given rate and CV for the QIF 

We will consider as given that the CV(1,D) {CV{—1,D)) is a monotonically increasing (decreasing) 
function of D; this wa s previously demonstrated by limit cases and by numerical evaluation of the integrals 
(jLindner et al.l . 120031 ) . Erom these properties we can conclude: for negative (positive) /i, decreasing 

(increasing) the CV from 1 (0) to 3"^^^, the parameter D changes monotonically from to infinity 
implying that each CV between and 1 has one unique contour line parametrized by the sign of /i and 
the value of D (see Eq.(|311)). The value ^ = is a special case where the CV attains exactly the value at 

the boundary between the regimes ^ < and A* > 0, namely, CV=3~^/^ (jSigeti and Horsthemkd . [19891 ) . 

If we can show that the rate or equivalently the mean ISI changes monotonically along the contour 
lines of the CV, then there is at most one intersection for a given pair of CV and rate and thus the 
mapping of rate and CV to and D is unique. Note that although our argument is similar to the one 
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Figure 5: Rate (a) and CV (b) as a function of the parameters /i and D for the QIF model. In (c) 
we show contour Unes for three values of the rate (dashed lines) and four values of the CV (solid lines 
with symbols) as indicated. The rectangular region in (c) is magnified in (d). In (d) we also plot the 
one curve (brown thick line) governing all ISI statistics in the Poissonian regime with Am(/j,, D) = 2.84 
(corresponding to the value &t D = 0.1,// = —0.36) resulting in /i = — (CTS-DAm)^/-^ (cf. the discussion 
around Eq. ((44)) ). The more sparse grid in the region close to /i = in (a, b) is due to the way we 
generated the points and does not reflect any property of the surfaces r{D, ^) and CV{D, 



used for the LIP, we will consider the change in the mean ISI along the curve of constant CV and not 
the change in CV (or variance) along a curve of constant rate as we did for the LIP. 
The directional derivative of the mean ISI along the CV contour line reads 
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From Eq. (|^ . we obtain: 



dD 



Ji^n-l/3n-2/3 

Im|3 



Using this expression and Eqs. (j37p and (|39p . one can rewrite the directional derivative as follows: 

2/3' 



V(T) • Vt 



1 2 

3L>^ ^ 1)2 da 



(47) 



(48) 



where we have used the auxiliary function I (a) from E q. (I37[) . The expression in the brackets vanishes 
by definition since it equals fi — ^cv according to Eq. (|42p . Hence we find that the last term is truly 
zero and thus the directional derivative of the mean ISI along the contour lines of the CV is negative 
throughout the {fi, D) plane 



V(T) - Vt 



1 



{T) < 0. 



(49) 



We have thus shown that (i) for each CV between and 1 there exists exactly one contour line and (ii) 
the mean ISI decreases always as we go along these contour lines in direction of increasing noise intensity. 
Hence, each mean ISI is at most represented once on a contour line of the CV and thus for each pair of 
rate and CV values there is at most one pair (/i, D). 



6 Conclusions 

To summarize, we have reviewed the behavior of rate and CV as functions of the input parameters for 
three different IF models. As the central result of our paper, we have shown that these statistics uniquely 
determine the input parameters for the models studied. This sets a framework for systematic comparison 
of these models: they can be compared on an equal footing when their parameters are tuned so as to 
yield the same rate and CV. Reports on these comparisons will be published elsewhere. 

It is tempting to consider the general IF model with white noise input: do rate and CV determine the 
input parameters for an arbitrary nonlinear function /n,odoi(i^) or equivalently for an arbitrary nonlinear 
potential U{v)7 Unfortunately, so far no general procedure to show the uniqueness of input parameters is 
known. What we showed in this paper relied on model-specific properties of the ISI moments for the PIF, 
LIE, and QIF. Approaches to the uniqu ene ss problem based on the general formulas for the moments of 
the first-passage time Eq. ([T5|l and Eq. could lead to conditions on the potential U{v) and the reset 
and threshold values of the IF model; however, we have not made any progress in this direction yet. We 
note that it may still be worth the effort to prove the uniqueness of parameter values for a given rate 
and CV for other specific neuron models of the IF type. One such a case is the exponential IF model 
( Fourcaud-Trocmc ct al., 2003), which has been successfully used to describe pyramidal neurons activity 
( Badel et all 120081 ). 

Another and more complicated open problem is to check whether the uniqueness of parameters also 
holds for more complex mode ls , i.e. multidimensional exte nsions of the simple IF mod el as those tak- 
ing ada ptation ILiu and Wand (120011) and thresho ld fatigue IChacron et al.l (120001. 12004 1) ; ILi ndner et al.l 

( 2005af) ■ subthreshold oscillations llzhikevichi (|2001l ) . or relative refractory effects lTuckwe ll ( 1978): lLindner and LongtinI 
( 20051 ) into account. In some of these cases, approximations to the firing statistics in certain regimes are 
known and so one can ask whether the additional variables change the uniqueness of parameters for a 
given rate and CV. 

Finally, the question of uniqueness of parameters could be also addressed in models that include a 
more realistic input beyond the diffusion approximation. This requires taking into consideration the 
shot-noise character of the synaptic i nput, i .e. th e fact that the neuron is driven by spike trains, as 
incorporated in Stein's model (Stcii:j, 119651 . 119671 ). The questions treated here may be worth to be 
addressed also for neurons which are subject to colored noise, either caused by a finit e synapt i c tirn e 
constant (jBrunel and Sergl 119981 ) or by temporal correlations in the pre-synaptic input (jLindneii . 120041 ) . 
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A Monotonicity of certain functions of interest for the LIF 

Here we want to prove that the two functions differences of which occur in Eq. (|36p arc monotonically 
decreasing functions of their arguments. Specifically, we argue that 

e°'erfc(a) > e'''erfc(6) if a < 5 (50) 

and 

erfc(a) Ji, erfc(o) 

We start by proving Eq. (|50p . For this purpose, we show that the difference of the left hand and the 
right-hand sides is always positive. Writing explicitely erfc(a) = dx exp(— x^) and performing 

the changes of variables s' = t — a and s" = t — h, we obtain for the difference 

\Jo Jo J V^r Jo 

if a<b, (52) 

since s > 0. Therefore, Eq. (|50() is proven. 

The second inequality Eq. (|5ip is equivalent to proving that the function 

= ° erfc(.) 

is a monotonically decreasing function of a. For this to hold true, the derivative of /(a) should be negative 
for all a, i.e. 

= ^ r dx ( e"'-'^'erfc'(a;) - e'''-"'erfc'(a)') < 0. (53) 

da 0Ferfc (a) Ja ^ ' 

The prefactor is positive and the integrand is strictly negative for all x > a. The latter can be seen by 

multiplying the i nteg rand by from which we obtain (e^ erfc(a;))^ — (e° erfc(a))^, which is negative 

by virtue of Eq. ([5Df . The proof of Eq. ([?T|) is therefore completed. 

B Some properties of the rate in the LIF model 

Here we show that for the white-noise driven LIF model the derivative of the rate with respect to the 
mean input is always positive. Further we calculate the limits of the rate for the mean input approaching 
minus and plus infinity. 

We first want to prove that 

drjdy. > 0. (54) 

Eqs. ^ and ^ imply that 

^{T)/^^i=-^{^{T)/^a + ^{T)/^b). (55) 



Now, 

We thus obtain: 



since a < b and by virtue of Eq. (|52|) . Eq. (|54p is therefore proved. 
Now let us prove that for the LIF one has 



d{T) /da + d{T) /db = ^{e^' erfc(6) - crfc(a)). (56) 
^^^^ - ' " (e"- erfc(6) - e"" erfc(a)) < 0, (57) 



dpi V 21? 




lim r = lim -— - = oo and lim r ~ lim -— - = 0. (58) 
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for which it suffices to show that the mean interval approaches zero or infinity as fj, goes to plus or 
minus infinity, respectively. The integrand in the integral expression for the mean interval Eq. (|24p is a 
monotonically decreasing function as shown above in Eq. (j52p : with this property we can estimate 

^{b - a)erfc(6)e^' < (T) < ^{b - a)erfc(a)e''' , (59) 

which is equivalent to 

erfc(6)e''' < \ < crfc(a)e'' . (60) 

V TT Vth - Vr 

Fo r /i ^ oo both a,b oo a nd th e functions on the left and right hand sides go to zero (sec eq.(7.1.23) 
bv lAbramowitz and Stegunl (jl970l) ). Thus, we obtain in this limit what proves the first of our limit cases 
in Eq. 

lim (T) = 0. (61) 

In the opposite limit of /i — oo, both a,b ^ — oo; the complementary error function attains a finite 
value in this limit (lim2;^_oo erfc(a;) = 2) and the exponential functions then yield a divergence of both 
sides yielding 

lim (T) ^oo. (62) 

fi — > — oo 

which proves the second of the asserted limit cases in Eq. ([55)) . 
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